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Abstract 

We present an event-by-event analysis of the cluster structure of final multi- 
hadron states resulting from heavy ion collisions. A comparison of experimental 
data with the states obtained from Monte Carlo generators is shown. The analysis 
of the first available experimental events suggests that the method is suitable for 
selecting some different types of events. 
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1 Introduction 



In many scenarios of heavy ion collisions, such as Quark Gluon Plasma (QGP)[1] or 
Disoriented Chiral Condensate (DCC) [2] formation, one expects an appearance of groups 
of particles, which differ somehow from the majority of produced hadrons. In particular, 
such groups are often characterized by small momenta in the rest frame of the group 
forming a sort of a cluster. Therefore it seems to be interesting to develop procedures 
looking for clusters in single events. 

We analyze the cluster structure of events resulting from the heavy ion collisions and 
from the Monte Carlo (MC) generators. In particular, we investigate, for single events, 
the dependence of the multiplicity distributions of charged hadrons in clusters on the 
parameter defining a cluster size. The procedure is presented in next section. We employ 
the cluster definition used in the procedure implementing the Bose-Einstein (BE) effect 
in MC [3]. We consider only clusters containing at least two hadrons and investigate 
the average multiplicity and the second normalized factorial moment, defining them with 
"two subtracted", i.e. n — 2, and 



This choice of moments is motivated by three facts. First, in our procedure each 
cluster contains by definition at least one hadron. Second, the percentage of "one hadron 
clusters" is not expected to obey any smooth multiplicity distribution (e.g. due to the 
contributions of "leading particles"). Therefore only the distribution for n > 2 may be 
expected to follow from some smooth formula. Moreover, if we consider the limit for 
F 2 at n — 2 — > for Negative Binomial Distribution in n — 2 with different values of 
the parameter k, we find it is finite and depends on k. Thus, in principle, our choice 
of moments allows to investigate the cluster structure in the limit of very small clusters, 
where other definitions lead to divergencies or to results independent of the shape of the 
distribution. 

In this note we work in a two-dimensional momentum space (rj — (j)) to be able to use 
the available experimental data of the KLM collaboration [4]. In principle, however, the 
data with full measurement of momenta are preferred, as they should allow to identify 
clusters better suited for the analysis of QGP or DCC signals. 

In the next section we present in more detail our procedure and the form in which 
the data are presented. We discuss also some modifications and generalizations of our 
approach. In the third section we compare the results of our analysis for the real - and 
MC generated events. The last section contains conclusions and outlook. 



To start the clustering procedure we have to define the distance in the momentum space. 
Since in the data set [4] only the pseudorapidity rj and the azimuthal angle are given, 
we use 



F 2 = 




(1) 



(n-2) 2 



2 Clustering procedure 



d 2 = (Arjf + (A0) 2 



(2) 
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as a distance measure. An introduction of different coefficients in front of the two terms 
in this definition (with their ratio between 0.5 and 2.0) does not change significantly the 
final results. Obviously, for data with full momenta measurement we would use some 
other measures, e.g. three momentum or four momentum difference squared. 

In our procedure originally each particle is considered as a single cluster. We fix the 
value of a "cluster size parameter" e and perform the clustering procedure for all pairs 
of particles. The particle is assigned to a given cluster if its distance d 2 to at least one 
particle from this cluster is smaller than e. Obviously for very small e values we have 
very few (or no) clusters with at least two particles. On the other hand, for large e values 
almost all particles must fall into a single large cluster. 

Our procedure is different from the clustering procedure applied in Ref. [4]. We do not 
define the "global cluster size" in the (77 — 0) space but consider only distances between 
particles. Moreover, the results are analyzed differently. In Ref. [4] one asks what is the 
percentage of particles in clusters; we investigate the multiplicity distribution of particles 
in different clusters for single events. The results are significant if we have a sufficiently 
large number of clusters (with at least two particles). This limits the possible range of 
e values, which must be neither too small nor too large. Thus our analysis is in fact 
complementary to the " intermittency analysis" performed in Ref. [4], where only very 
small bins in momentum space are relevant (corresponding to small e values). At the 
lower end of our range most of the clusters contain just two particles; at the higher one, 
as we shall see, there are typically a few tens of particles in a cluster. To characterize the 
muliplicity distribution we calculate for each event the average multiplicity of a cluster 
and the second scaled factorial moment, as defined in the Introduction. 

3 Results: comparison of MC with data 

We performed the clustering procedure for the four KLM events presented in Ref. [4] and 
for the events obtained from the VENUS generator [5]. In addition we use the "random 
events" of similar multiplicities obtained with a uniform random generation of (77, 0) points 
(SERENE events). The results are shown for the range of e values for which there are at 
least 40 clusters (with more than one particle) in the event. We checked that the results 
do not depend on the order in which the data are read in, as expected. 

The average cluster multiplicity (or, to be more precise, its excess over the minimal 
value of 2) increases monotonically with e. It is interesting to note that the value of e 
for which the increase becomes faster than linear is the same as a value for which the 
number of clusters with at least two particles starts to decrease (for smaller e the decrease 
resulting from joining small clusters into the large ones is compensated by formation 
of small clusters from single particles). As shown in Fig. fa, the increase is faster for 
events with larger global multiplicity N. Two lower data sets correspond to N — 742 
and N = 743, and the two higher ones to N = 926 and N = 986, respectively. If the 
distance e is scaled by the average distance e between neighbouring particles (inversely 
proportional to N, e ^ 12tt/N), differences between events disappear as shown in Fig.fb. 
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Fig.l. Tie average hadron multiplicity of a cluster n — 2 (Tor clusters with at least 
two particles) as a function of a)the parameter e defining the maximal distance in the 
(rj — (f>) space, for which two particles are joined into a cluster, b)the same parameter 
scaled by e = 12n/N. Data are taken from Ref. [4]. Four kinds of symbols correspond 
to four different events. 

In the MC generated events a similar increase occurs. The corresponding plots (not 
presented here) show no clear differences between the experimental and VENUS events, 
although the latter ones have typically higher global multiplicities resulting in higher 
average multiplicities of clusters. In random events there seems to be a larger spread of 
data points for fixed global multiplicity, but again no clear difference shows up. 
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Fig. 2. The second scaled factorial moment (Eq. (1)) as a function of average multi- 
plicity n — 2 for the events obtained from the VENUS generator. Each symbol corresponds 
to one event. 

On the other hand the analysis of the second scaled factorial moment as a function of 
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average multiplicity shows a striking differentiation of events. For small n — 2 (below 4) 
it rises irregularly to the values of a few units. Then in all the VENUS events we see a 
much faster rise to the values of about 30 - 40 as shown in Fig. 2. For transparency only 
five out of twenty analyzed events are shown. 




Fig. 3. The second scaled factorial moment as a function of average multiplicity n — 2 
for the events obtained from a uniform random numbers generator. Each symbol corre- 
sponds to one event. 

In the random events the rise is much slower (Fig.3), although we selected multiplicities 
corresponding to the maximal multiplicity of the VENUS events. Out of four KLM events 
two seem to follow the "VENUS pattern" and two look like "SERENE events", as seen 
in Fig.4. 
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Fig.4. The second scaled factorial moment as a function of average multiplicity n — 2 
for the KLM events. Each symbol corresponds to one event. 

The irregularity of F 2 for small n — 2 shows that for the events with global multiplicity 
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of the order of 10 3 we cannot draw any significant conclusions about the limit n — 2 — > 0. 
Our choice of the definition of F2 "with 2 subtracted" (1) is found a posteriori not too 
relevant. However, the situation may be different at higher multiplicities expected at 
RH1C energies. 

One should stress that large values of F 2 signal large fluctuations of the cluster mul- 
tiplicities (as compared to the average value). The violent growth observed in VENUS 
events suggests an occurence of one or a few clusters with multiplicity of the order of 100 
when the average fluctuates around 10. Apparently this does not occur in the randomly 
sampled events and only two out of four KLM events show such behaviour. It is interest- 
ing to note that these two events are the same which show fast " intermittent" rise of the 
scaled factorial moments in small bins of the phase-space; this rise was interpreted as a 
manifestation of "strong dynamical correlations" [4]. 

It is, in fact, rather surprising that the same events seem to be "special" in the limit 
of very small bins (" intermittency analysis") and for a quite large average cluster size in 
our analysis. 

Even more surprising is the fact that in our analysis all the "VENUS events" look like 
"special" KLM events. It may be so that fixing the parameters of the VENUS generator 
by fitting the averages over many experimental events induces such cluster structure in 
all single events, whereas the data are more differentiated. It will be really interesting 
to repeat this analysis on a much larger sample of events (both experimental and MC 
generated) and for other choices of variables and clustering procedures. 

4 Conclusions and outlook 

We have applied a simple cluster analysis to single events of heavy ion collisions, both 
experimentally and MC generated. Analyzing the multiplicity distributions of charged 
hadrons in clusters we found two different patterns of behaviour already in the first four 
published experimental events. On the contrary, all the VENUS generated events follow 
one pattern. Thus the reasonable agreement of MC description with the inclusive features 
of data (obtained by averaging over many events) does not apply for single events. 

We conclude that the cluster analysis of single events may allow to select special classes 
of events signaling, e.g., the formation of some exotic states (QGP, DCC). Obviously our 
results present only a preliminary evidence; similar analysis should be performed using 
different variables. In particular, data with full momentum measurements and particle 
identification should be compared with the MC results. Different clustering procedures 
should also be introduced. The results may help to improve our understanding of multiple 
production in heavy ion collisions. 
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